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We compute the second-order matching conditions for tensor metric perturbations at an abrupt 
change in the equation of state. For adiabatic perturbations on large scales the matching hyper- 
surface coincides with a uniform-density hypersurface. We show that in the uniform- density gauge 
both the tensor perturbation and its time-derivative are continuous in this case. For non-adiabatic 
perturbations, the matching hypersurface need not coincide with a uniform-density hypersurface 
and the tensor perturbation in the uniform-density gauge may be discontinuous. However, we show 
that in the Poisson gauge both the tensor perturbation and its time-derivative are continuous for 
^vj adiabatic or non-adiabatic perturbations. As an application we solve the evolution equation for 

second-order tensor perturbations on large scales for a constant equation of state and we use the 
matching conditions to evolve the solutions through the transition from an inflationary era to a 
radiation era. We show that in the radiation era the resulting free part of the large-scale tensor per- 
turbation (constant mode) is slow-roll suppressed in both the uniform-density and Poisson gauges. 
Thus, we conclude that second-order gravitational waves from slow-roll inflation are suppressed. 

PACS numbers: 98.80.-k, 04.30.-w, 98.80.Cq 

U 

Recent precise measurements of the Cosmic Microwave Background (CMB) anisotropies have revealed the nature 
. of the primordial perturbations (l| and future experiments such as Planck [2] will improve these measurements 
significantly. However, CMB anisotropies and large scale structure probe density perturbations only on large scales 
and there is no direct way to observe density perturbations on smaller scales except for the case where large density 
l— "™ 1 ' perturbations form primordial black holes [3, H| • Recently, it has been recognized that gravitational waves can be 
used to probe small scales physics like preheating 0-0| and density perturbations on small scales [1, At first 
order, tensor perturbations decouple from density perturbations. However, at second order, density perturbations 
can generate tensor perturbations (loj . The second-order tensor perturbations have been calculated during inflation 
, [ll| in the radiation era [8] and in the matter era H3 - ll7j . While the amplitude of the second-order tensor 
perturbations is generally very small (determined by the square of the amplitude of density perturbations) , there are 
several interesting situations where the second-order tensor perturbations give interesting observational consequences, 
t^- ■ For example, if density perturbations are large enough to form primordial black holes, then these density perturbations 
can generate sizable second-order gravitational waves @. In the matter era, the second-order tensor perturbation 
remains constant, which enhances the power spectrum of the tensor perturbations This happens either in the 
late-time matter era or in an early-time matter era due to the oscillations of the inflaton field after inflation [l7^ . 
The first-order tensor perturbations generated during inflation are model dependent and their amplitude can be very 
small. However, the amplitude of the second-order tensor perturbations is determined completely by that of the 
density perturbations. Thus this gives a lower bound for the tensor to scalar ratio [H[. 

There are still several important issues in the second order generation of gravitational waves. At second-order, 
the tensor perturbations are not gauge invariant under the first-order gauge transformations [l9| . Thus they depend 
on the choice of the gauge for first-order perturbations. It is however possible to construct gauge-invariant tensor 
perturbations by eliminating the gauge degrees of freedom [20H241 ] . A commonly used gauge is the Poisson gauge 
[HI EH, [25| and almost all calculations of second-order gravitational waves have been done in the Poisson gauge. 
During inflation, another commonly used gauge is the uniform-density gauge. This gauge is often used to calculate 
higher-order correlation functions for the curvature perturbation (2o| and these calculations involve the second-order 



tensor perturbations because they are coupled to the curvature perturbations at higher order [27J. In fact, the 



amplitudes of tensor perturbations in these two gauges can be very different. In the Poisson gauge, the amplitude 
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of the second-order tensor perturbations during inflation is suppressed by slow-roll parameters but it is not in the 
uniform density gauge. 

The question is then whether the amplitude of gravitational waves after inflation is slow-roll suppressed, or indeed 
whether large tensor perturbations can be generated at second order by the transition to a radiation dominated 
universe. In order to address this problem, we should track the evolution of the tensor perturbations during the 
transition from inflation to the radiation era. 

At first-order, matching conditions for cosmological perturbations have been developed in cases where there is 
an abrupt change in the expansion rate at the matching surface |28l - [3l| . such as reheating at the end of inflation. 
In this paper, we extend this matching condition to second-order perturbations and derive matching conditions for 
second-order tensor perturbations on large scales. Matching conditions will be derived in two cases. For adiabatic 
perturbations on large scales a sudden transition in the equation of state must occur at a specific density. Thus 
we apply the matching conditions across a uniform-density hypersurface. We will also consider the case where the 
transition surface is given by a fixed value of some field x, which need no longer coincide with a uniform-density 
hypersurface in the presence of non-adiabatic perturbations. In this case the matching conditions are applied across 
a uniform- x hypersurface. These matching conditions are applied to the second-order tensor perturbations generated 
during inflation and we calculate the second-order tensor perturbations in the radiation era after reheating in both the 
Poisson gauge and the uniform-density gauge. We will show that the constant mode of the second-order gravitational 
waves on large scales, corresponding to the free part of the tensor perturbations, is the same in the two gauges in the 
radiation era and it is suppressed by slow-roll parameters. 

This paper is organized as follows. In section II we summarize the gauge transformations at the first and second- 
order and derive solutions for the metric perturbations on large scales in the Poisson gauge and in the uniform-density 
gauge. In section III, the evolution equation for the second-order tensor perturbations is derived. General solutions for 
the second-order tensor perturbations are derived on large scales with a constant equation of state and we check the 
consistency of solutions in two different gauges under the gauge transformation. In section IV, matching conditions 
for perturbations are developed at the first and second-order in the cases where there is an abrupt change in the 
expansion rate at the matching surface. We consider both adiabatic and non-adiabatic matchings. In section V, 
we apply the matching conditions to the second-order tensor perturbations generated during slow-roll inflation and 
calculate the evolution of the second-order tensor perturbations on large scales in the radiation era. Section V is 
devoted to conclusions. 



II. GAUGE-INVARIANT PERTURBATIONS 



In this section, after introducing the necessary notation and the perturbed metric, we shall review the first and 
second-order general gauge transformations for the perturbations. After that, using some simplifying assumptions, 
we will obtain the gauge transformation rules for second-order tensor perturbations. We will also derive the solutions 
for first-order scalar metric perturbations on large scales. 



A. Perturbed metric 

Throughout this work Greek indices /i, v, . . . take values from to 3, while Latin indices denote spatial 

indices and can take values from 1 to 3. The components of a perturbed spatially flat Friedman- Robertson- Walker 
metric (FRW) can be written as 



5oo 



9oi 



9ij 



= - a (,) 2 (l + 2gi^), 

C+oo \ +oo 1 

r=l ' / r=l 



(1) 

(2) 
(3) 



where rj denotes conformal time and Cy is traceless, i.e. C£ = 0, where the index was raised with <S* J '. The 
superscript (r) indicates the order of the perturbation. The perturbations depend on conformal time and position x. 
The inverse metric up to second order can be found in Appendix. We decompose the perturbations in scalar, vector 
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and tensor parts as 

Bl r) = diBV+S$ r \ (4) 



where d'Sf ' = 0. 

C^=2D ij E^+d i F^+d j Ft ) + h ( ^, (5) 

where d % F^ r ' = and h 1 ^ — d l h[^ = 0. The operator is defined as = didj — ^Sijd 2 . The energy density p 
and the four velocity of matter are perturbed like 



p=p(°)+J2W r) > ( 6 ) 



r 

r=l 



(7) 



u M is normalized as u M u M = — 1 and this implies that related to the lapse perturbation A^ r \ The velocity 

perturbation is decomposed as 

v i{r) _ d i v \\{r) +u Xi(r) (g-j 

B. Gauge transformations 

Bruni e£ aZ. [19( have shown us that given a tensor T, its perturbations in two different gauges X, and Y are related 
at first order as 

5T Y{1) -6T x(1) = £f ) T(°\ (9) 

and at second order we have 

(£^> + £ 2 ei) )T^ +2£ iW ST x{1 \ (10) 



ST Y(2) - ST x{2) 



where £W and £^ 2 ) are the gauge shifts at first and second order respectively. The Lie derivative along a vector field 
£ of a rank two covariant tensor T is given by £{T^ U — T^ Vt \^ x + £^7a„ + T m a- We also decompose £^( r ) is scalar 
and vector parts as 

fO(r) _ a (r) ) £t(r) = 5 i^(r) + d t(r) ) 

where did 1 ^ = 0. 

J. First order 

At first order, the gauge transformations, Eq. ([9]), can be written as 

I« = + (12) 

C« = C^-i^W-WaW, (13) 

^ = itf> -««+/#>' + «?>', (14) 
C« = cW + ZD^W+da+dS, (15) 
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= B« 






= s« - 






= 




i 


i 




h {1) 







where prime denotes the derivative with respect to conformal time and the conformal Hubble rate is % = — . The 
last two equations imply 

(16) 
(17) 
(18) 
(19) 
(20) 

Four dimensional scalars, like the energy density or the pressure, have the transformation rule 

5p = Sp + p (0) 'a (1) , SP = SP + P^'a^. (21) 
The spatial component of the four velocity transform as 

„m = v iW- d ipm' - d m', (22) 

or equivalently 

i m =v KD_pW' ) c xi(i) =v x i (i)_ d i(i)' j (23) 

while the temporal component transforms as because of the normalization constraint ti 0( 4) = — A^K 

Following Bardeen [32|], we can define gauge invariant quantities by specifying completely the choice of coordinates. 
The longitudinal or Poisson gauge is defined at first order by requiring Bp^ = Ep = 0, which from Eqs. (jTHJ) and (Ti~8")) 

fixes Pp" 1 — —E^ and = — E^', where the subscript P denotes the Poisson gauge quantity. The remaining 
scalar metric perturbations then have the gauge invariant definitions 

$ ee A$ =A& +h(bW -E^') + (b^ -E W )' , (24) 



* ee C P 1] = C« - H (b™ - + l -d 2 E^\ (25) 

The gauge- invariant density perturbation corresponding to density perturbations in this specific gauge is, from Eg. (1211) . 

Spf =5 P + P ^'(b^ -E^') . (26) 



An alternative set of gauge invariant variables can be constructed by working on uniform-density hypersurfaces for 

(i) 

UD 



which SpuD = 0, which from Eq. (|2"Tj) fixes ajj' D — —Sp/p^'. The curvature perturbation in the uniform density 
gauge is thus the gauge invariant variable 



- C = C« + \d 2 ^ ] D = CW + \d 2 E^ + -J^p« (27) 

Here and in the rest of this work we have choose the threading of the time slices, /3^ D = —E^\ such that Ejp D = 0. 
The gauge invariant definition of the scalar part of the shift vector (|16p in the uniform density gauge is given as 



p(oy p (o)/ ^ n 



In the absence of anisotropic stress, the Einstein equations in the Poisson gauge require $ = \& and the energy 
constraint equation gives a relation between the curvature perturbation in the uniform-density and Poisson gauges 

r=- m ~ H - + - d 2 ^ (29) 

^ n 2 -w u 2 -w ^{n 2 -w) ' y ' 

In the next section, when writing the equations of motion for second order tensor perturbation, we will neglect terms 
with more than two spatial gradients. We also ignore non-adiabatic perturbations. This implies that we can take 
C = at zeroth order in the gradient expansion. The scalar part of the Einstein equations gives 

a- 2 (6H(C {iy + HA<») - 25 2 C«) = -K 2 G 6p<», (30) 
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where kg is the gravitational constant. In the uniform-density gauge Sp^ = 0, thus on large scales Aue, — 0. For 
a perfect fluid with a constant equation of state w we have £ = V =constant on large scales, and Eq. (j2U]) gives the 
solution 

9 = -^V + Wv-^, (3D 
5 + 3w 

where W = W(x) is a constant of integration to be determined by the first order junction conditions. If w — 1/3, the 
previous equation simplifies to give 

= -\v + Wr]~ 3 - (32) 



2. Second order 



The second-order gauge transformation can be found in Bruni et al. For the trace-free spatial metric perturbation 
it reads 



(7 (2) = c (2) 



+2 (-4C« + a^do + e W d k + fflat") (djg) + 
) 



+2 
+2 

+ 2 (d^j+Dy^W 

We can define gauge invariant second order tensor perturbations in the Poisson gauge as 

h$> P = ^ ) +2[(c«' + 2^W)aW + cgU fe(1) 

+ (-4C« + a^d + e {1) d k + maW) (dg } + DyfiW) 

+(<-«S ) +€§ ) >S+Kii+^+4lW€a ) 



(33) 



TT 



(34) 



where all the quantities in the r.h.s. of the previous equation are in an arbitrary gauge and the parameters of the 
gauge transformations are explicitly given as — — = —E^\ d^' — —F^ and = d l j3^ + d 1 ^ . 

In a similar way, a gauge invariant definition of the uniform density second order tensor perturbation is 



■ (2) 

ij UD 



TT 



(35) 



where a' 1 ' = -SpW/pW, = df ] = -F z {1) and = + <* i(1) . 

In this work, we will neglect the effect of first-order vector and tensor perturbations in order to simplify the second- 
order equations as much as possible. In many models of the very early universe, such as slow-roll inflation which 
we consider later, vector and tensor perturbations are suppressed with respect to scalar perturbations. First-order 
primordial vector and tensor perturbations, if they do exist, would provide an additional, independent source of 
second-order gravitational waves. 

We choose to set E^ to zero in both Poisson and uniform density gauges. This implies that in the gauge transfor- 
mations between these two gauges the parameter is zero (and ef ^ = 0). The first order temporal gauge shift, a*- 1 -*, 
is given in Table I. For example, if we start in the Poisson gauge and if we do a gauge change with = —Sp^ / p^ 
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TABLE I: Time gauge shift cr 1 ' at first order between our different gauges. 





Poisson 


Uniform density 


Poisson 







Uniform density 


B W 






then we obtain the quantity in the uniform density gauge (tilde gauge), as given in Eqs. (|12j|15[) . The Poisson gauge 
is defined by B^ = E^> = F^ = and the uniform density gauge is Sp^ — — — 0. With the previous 
assumptions and with the choice E^> = 0, the gauge transformation equation (|33p for the traceless perturbation C 



(2) 



simplifies considerably to 



C 



(2) 



a 



(2) 



2B 



(i) 



a 



J 1 ) 



,;) 



which gives the following rule for transforming second order tensors 



h {2) = 



,(2) 



,3) 



TT 



(36) 



(37) 



where TT means the transverse and traceless part of the terms in between square brackets. In the case of transforming 
from the Poisson gauge to the uniform density gauge the previous equation further simplifies to give 



h {2) 

a ij UD 



l {2) 
b i3 P 



2 



[(£ + *)(£ + *) 



(38) 



The right-hand-side of this equation leads to different solutions for second-order tensor perturbation in uniform-density 
or Poisson gauges. Nonetheless we have constructed gauge-invariant expressions, (|34l) and (|35|) . for the second-order 
tensor perturbation in the uniform-density and Poisson gauges, and the difference between them is also manifestly 
gauge-invariant . 



III. TENSOR PERTURBATION EQUATIONS 



In this section, we will solve the first and second-order large-scale equations of motion for tensor perturbations in 
a fluid with a constant equation of state, in both the Poisson and the uniform-density gauges. At the end of the 
section, we shall present the second-order gauge transformations and the relations between the integration constants 
that appear in the two gauges. 



A. First order 



The first order tensor perturbation is gauge invariant and its equation of motion on large scales is: 

h$" + 2Hhf = 0, (39) 

that can be easily solved 

hf) =Y ij +Z ijV -^, (40) 

where and Zij are transverse and traceless integration constants to be determined by the initial conditions. During 
matter or radiation eras, Z+j is a decaying mode and the growing mode is constant in time, as expected for first-order 
gauge invariant tensor perturbations on super-horizon scales. 



B. Second order 



In this subsection we will present the equation of motion for second order tensor perturbations in the uniform 
density gauge and in the Poisson gauge. Note that we neglect first order vector and tensor perturbations and we 
choose E^ = 0. We also ignore first and second order anisotropic stress. 
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Using the second order Einstein equations from Bartolo et al. review [33| we get 



where 



T 



i(2 



TT 



(p<°> +PW) [c^W +djB^) 



TT 



(41) 



(42) 



G 



i(2 



TT 



-2 



d i AdjA + 2Ad i d j A - 2Cd l d 3 A - - fl^C + S&CdjG + 4C<9 i <9 i C 



+2rld l Bd 3 A + ArlAd l d 3 B + A'd l d 3 B + 2Ad i d j B' + d 2 Bd i d j B - d 3 d k Bd i d k B 
-ZH&CdjB - 2Hd l BdjC ~ d l C'd 3 B + d 3 C'd l B - d l Cd 3 B' - d 3 Cd l B' - 2Cd l d 3 B' 

n TT 

+C'd l d,B - mCdi&B 



(43) 



where in the previous equation we omitted the superscript (1) to indicate the order of the perturbations and we 
corrected the typos of equation (A. 43) of Ref. (33[. 

1. Poisson gauge 

The equation of motion for second-order tensor perturbations in the Poisson gauge is known and can be found in 
Baumann et al. [HI (see also [1, [H, [34| ) . Under our assumptions (for constant w) it simplifies to 



T di (4-' + W$>) d 3 (*' + H9) 



TT 



3(l + u;) H 2 

Using the solution for ^ as in equation pip, the equation of motion can be written in a simple form as 



,(2)' 



= 8 



3 + 3w 



[d i Vd j V) TT + 4(5 + 3w) [d.WdjW] 1 1 t]' 2 ^ 



r P r P n 5 + 3n 



(44) 



(45) 



where we drop the gradient terms on the left- hand-side of (1441) . This is the equation of motion for second-order tensor 
perturbations, at leading order in a gradient expansion. It can be solved to give 



fctfP = YiiP + Zijpr)-^ + 4 ^ T k 7Za2^' [diVdjV]" rf + [BiWdjWY 1 rj-^, (46) 



(l + 3w)(3 + 3w) 



iTT 2 



(l + 3w) 2 



n TT 



(5 + 3w) 5 



where Yi 3P and Zi 3p are integration constants. 



2. Uniform density gauge 

Using Einstein equations (|41j) . we can derive the equation of motion for second-order tensor perturbations in the 
uniform-density gauge as 



i(2)' 



= -45 

J UD V 



TT 



(47) 



where the source is given by 



S l3 = d l Ad 3 A + 2Ad i d j A + 2(d i d j A + d 3 Ad\ + d l Ad 3 ( + 30*C0jC + AC&djC + 2Ud i Bd j A + AUAd l d 3 B 
+A'd t d j B + 2Ad l d 3 B' + d 2 Bd l d 3 B - d 3 d k Bd l d k B + 2UdX'd j B + 2Hd l Bd 3 ( + d\'d 3 B - d 3 ('d l B 



+d l C,d 3 B' + d 3 (d l B' + 2(d l d 3 B' - £d l d 3 B + mC,d 3 d l B - -fy (-(' + HA) cT 2 <9 ! ( C' + ^ - — £P 



p(o) + p(0) 



(48) 
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and we have used the first-order Einstein equations to simplify the result. Using Eqs. (|28|) . (|3Tj) . Aud — on 
large scales, ignoring higher derivatives terms and the gradient term on the left-hand-side of (|47l) . we can write the 
equation of motion (|47|) for second-order tensor perturbations in the uniform-density gauge at leading order in a 
gradient expansion as 



$ ) l D +^ ) ' UD =mVB j Vf T . (49) 



This can be easily solved 



h f UD = Y ^ud + ZiiunV- 3 ^ + 2^|^ [tkVdjVT?, (50) 



where Yij UD and Zij UD are integration constants to be determined by the initial conditions. 
The second-order gauge transformation ([38)1 between these two gauges is 

,(2) ,(2) ofl±^Vr*T/*.T/lTT 2 (I+JH 2 r 0ma ,^TT. 



K: ? UD = W P -2 [^-^ ) [diVd jV y ^-L- [WW f,-r** 

[d . wd v + d 3 w d%V f T (51) 

5 + 6W 

We have checked that Eqs. (|49|) and (|45j) are related by (|5"Tj) (ignoring higher gradients). This provides a good 
consistency check on our calculation. In a similar way, one can confirm that the solutions of Eqs. (H9l and (l4"5l) arc 
related by ([5"Tj) . and we get two relations between the integration constants in the two gauges as 

Y — Y 

Z lJp = Z lJUD + { \ + 3w)2 [diWdjV + djWdiV] TT . (52) 



We note that the gauge- dependence of the second-order tensor perturbation, (fBlj) . affects only the amplitude of the 
time-dependent parts of the perturbation in Poisson or uniform-density gauges, and leaves the growing mode, Yy, 
unaffected on large scales. 

IV. MATCHING CONDITIONS 

In this section, we shall derive the matching conditions for the metric and the extrinsic curvature up to second 
order. These matching conditions are applicable on large scales and when there is an abrupt change in the expansion 
rate at the matching surface. We will first consider the case where this transition happens at a specific energy density, 
as it does for adiabatic perturbations. We will then discuss the more general situation, where the transition is not 
determined by the density but by some other field X- We call this case the non-adiabatic matching. (See section 9 of 
Ref.[24[ for the definition of non-adiabatic perturbations of multiple fluids or fields.) In the final subsection, we will 
apply these matching conditions to the second-order tensor perturbations coming from first-order scalar-perturbations. 

Matching the induced metric and extrinsic curvature of a uniform-density or uniform-^ hypersurface specifies a 
physical model for the matching, but does not specify the gauge. We will choose to work in a gauge in which the 
matching surface coincides with a constant-?? hypersurface. This is not a physical restriction; it is simply a choice of 
gauge, reflecting our freedom to choose the time-slicing. 

The metric of the constant rj hypersurfaces (in the perturbed spacetime) is 



Hj = 9ij = a{r)f 



r=l ' / r=l 



(r) 



(53) 



The unit time- like vector field orthogonal to these surfaces is iV M = (N°, N l ) with the normalization N^N^, = — 1. 
is given by q^ v — g^ + N^N^ which together with equation (I5"3l implies = (N , 0). The inverse induced metric 
is q»" = g^ v + N»N V and it gives 

q n(o) = g n(o) ) q im = g m) ) q im = gtm + a - 2 &v>B*v>. (54) 

The extrinsic curvature of the constant r\ hypersurfaces is 

= 2 £ N*<ln» = 2 (W-^ + K^™ + N U^) ■ ( 55 ) 
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Defining the trace and the traceless parts of the extrinsic curvature like 



K = q lJ K lh Kij = Kij - -qijK, 



(56) 



the junction conditions are 



[g«]+ = 0, [Kijt = 0, 



(57) 



which are equivalent to 



[g«]i = 0, Kij =0, [K}t = 0. 



(58) 



A. Background 



In the background, 



and 



Therefore 



N°M = T a-\ N^=0, N^=±a, (59) 



4°) = T a'%. (60) 



#(°)= T 3-, K^=0. (61) 



Continuity of the induced metric and its extrinsic curvature, Eq. (|57|). at zeroth-order thus requires that a and a' are 
continuous, i.e., 

[0)1 = 0, [H]t=0. (62) 

These matching conditions can be used to give the initial conditions for the scale factor and its derivative after a 
sudden transition in the equation of state from a u>_ era to a w + era (such as from inflation to a radiation domination 
era). The solution of the background Einstein equations for the conformal Hubble rate is 

n=( 1 -±^r 1 + C 1 Y\ (63) 



2 

where Ci is an integration constant. From the continuity of the scale factor and its first derivative we can obtain 

n = f r 1 + C l _\ for 77 < 77*, 

n = / l + 3w + ri + c \ for TfrKn, (64) 



where the constant C% . is 



We can write the previous equations as 



3 

C 1+ = C\_ + -?7* (w- - w + ) . (65) 



, T + 3w— \ 1 
" = I 2 ^"J V-<V-*, 

H = l+p±\ 1 for ^ < (66) 
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where the new time variables rj± are 



<1± = <1 - 



1 + 3w± 



Cl±! 



(67) 



and the transition times in the new variables are 

2 



V-* = V* 



-Ci_, 



1 + 3w. 



-V* 



-C x _. 



1 + ,T * l + 3w+" l + 3w 

The subscripts + (— ) denote that the quantity should be evaluated at a time after (before) the transition time rj*. 



(68) 



B. At first order 

The first-order matching conditions have been calculated previously in Refs. [29l - l3l"l ] . 

Continuity of the induced metric at first order leads to the conservation of the spatial metric perturbations CW 
and C[]\ or equivalently 



F, 



(i) 



= 0, 



0, 



E 



(i) 



= 0, 
= 0. 



The first-order correction to the ortho-normal vector field is 



and the extrinsic curvature is then 



K 



(i) 



± a - 
2 



26; 



(2UC^ + HA& + c W"j + B W + B W _ 2HC[f - 



or equivalently 



K 



(*>i) 3 



2 *■» 



(69) 

(70) 
(71) 

(72) 

(73) 
(74) 
(75) 
(76) 



Note that Eq. (|73|) combined with Eqs. (|69|) and the Einstein equations, Eq. (|30|) . enforces continuity of energy density 
across the hypersurface 



The junction condition (|57[) for the extrinsic curvature gives 







+ 


-o, 


5(1) 


_£(!)'' 


+ 


= 0, 




j 


+ 


= o, 






i + 


= 



Sp 



i:| + =0. 



(77) 



1. Adiabatic matching 

For adiabatic perturbations on large scales a sudden transition in the equation of state must occur at a specific 
density. Thus we apply the matching conditions across uniform-density hypersurface [29| . We can use the freedom in 
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the choice of spatial coordinates on the matching hypersurface to set = 0. In this case the matching conditions 
and (j74")l across a uniform-density hypersurface reduce to 

[C]+=0, [¥]t=0. (78) 

Note that the remaining junction condition for scalar perturbations, Eq. ([75)1 , reduces to = on large scales 

which is automatically satisfied for adiabatic perturbations since £' = 0. Taking the growing mode for adiabatic 
perturbations on large scales before the transition 

3 + 3w_ 

= V- = const , (79 

then the junction conditions (|78p. together with the general solution for constant w, Eq. (|3"Tj) . can be used to determine 
the solution of ^> after a sudden transition as 

3 + 3w4 



*+ = - , ' „ Z- V+ + W+Tt + ~ m + , (80) 
5 + dw + 



where V+ = V— = C- an d the amplitude of the decaying mode, W+, is given by 



= 3 + 3^ ^ / _ 3 + 3._5 + 3. + y . 
5 + 3w+ \ 5 + 3iu_ 3 + 3iu_|_ / 



3. N on- adiabatic matching 



To consider a more general matching condition we will consider the case where the transition surface is given by a 
fixed value of some field, x- The scalar metric perturbations on a uniform-^ hypersurface can be given in terms of 
the metric perturbations and x-field perturbations in an arbitrary gauge as 



(i) _#(!)' = r« _ pay i 



i^ 1 



/(0)' 



Departures from adiabaticity are characterised by perturbations of the %-field with respect to the total density 

V 1 ) 5 X (iy 



S X = H 



,(0)' 



X 



(0)/ 



(82) 
(83) 

(84) 



We will recover the adiabatic matching conditions when S x = and x is unperturbed on uniform-density hypersurfaces, 
so that the uniform-^ hypersurface coincides with a uniform-density hypersurface. 
The usual gauge-invariant curvature perturbations can then be given as 



C — Cx ~ $x > 



x x 

The junction conditions for scalar metric perturbations (|69l) and (IT4[) require 

B?) - si 1 )'] 



[C x ]I = o, 



0. 



(85) 
(86) 

(87) 



while Eq. (f73"| enforces energy conservation across the hypersurface, Eq. ([77)1 and thus 



.1 + w+ 

Expressing Eqs. (|87[) in terms of the usual gauge-invariant curvature perturbations we then have 

1 + W- 



[C]I= i- 



l + w+ 



S 



x- ' 



$(i) 



= 0. 



(89) 
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Analogously to the previous subsection, these junction conditions can be used to evolve the solution of "J/ across 
the transition. The solution is given by Eqs. (f?9"f and f|5U|) . where now the integration constants after the transition 
are given by 



1 + W- 
l + w+ 



(V--&-). W+ = 



3 + 3wj 



V + - 



3 + 3w- 



5 + 3w+ 5 + 3w- 



V- V 



5 + 3m + 



(90) 



C. At second order 

Continuity of the induced metric at second order leads to the obvious extension of the first-order junction conditions 
69| for the spatial metric perturbations C' 2 ' and , or equivalently 



C (2) 



(2) 



= 0, 



,(2) 



At second order, in the ortho-normal vector field we have 

-i f_ A (i) B m + 2C^B^ - B^CPW + 



N o(2) = ± ±_^ 
2 



N*P> = ± a - 



The extrinsic curvature tensor is 



0, 
= 0. 



-± a ~ 

2 

_B»(2) 



(91) 

(92) 
(93) 



A" 



(2) _ 



UC V) + -C^' + -UA^ - 2UA^C^ - A«C«' + -u( - 3 (a^Y) - B^C™ 



± a 



2 n 4 u 



(2) 



ij,k 



CM) >(» j) fe (j> ; 



(94) 



The trace and the traceless part read 

K w = ±JL [ 3 (n A (2) + C (2)') + B m + 2 ( 3C ,(i)' + B ,m^ ^ 2C w _ A ^ _ 9 ^ (^(d)' 
_ 2C ,,i(D o(d + mB v B m _ 2 C tf < 1 >B# - 2B*«qf } + c ij ^c[f 



K 



(2) 



±a 



C 



(2)' 



2 fort 2 



-AMb£\ + 2C^B^ - B^C™ - B^C^l + B«W}1 
(*>j) >(* J) M 1 J) fc (*>j) A 1 J) fc 



3 V 2 

The remaining matching conditions are then 



x / n fe(2) -. N 



'k 



'kl 



K 



(2) 



(95) 



(96) 
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D. Matching conditions for second-order tensors from first-order scalars 

1. Adiabatic matching 

For adiabatic perturbations we match on a uniform-density hypersurface, Sp = 0. Matching the induced metric on 
the transition surface then implies, from Eq. (|9ip . that 



,(2) 

y ud 



0. 



(97) 



Setting — and neglecting first-order vector and tensor perturbations, so that = 0, the matching condition 
for the transverse and traceless part of the extrinsic curvature (1^51) gives 



■ (2)' 

y ud 



i + 



= -4 



( A UD B UD,i 3 + K(i B U. 



UD,j) 



TT 



1 + 



(98) 



Continuity of C and B^Jj across the matching hypersurface is enforced by the first-order scalar matching conditions 
(IMl) and (j?4")) . In addition, in the uniform-density gauge, from ([73]) we have 



ri 



(99) 



Thus for adiabatic perturbations on large scales, for which (' = 0, we find 



(2)' 
ij UD 



= 0. 



(100) 



Thus the tensor metric perturbation and its first derivative are continuous in the uniform-density gauge for adiabatic 
perturbations on large scales. 

The second-order gauge transformation from the uniform-density gauge to the Poisson gauge yields 



h m = h {2) 

l 3 UD l 3 P 



2(5^ 



,(2)/ 



k ud y p 
For adiabatic perturbations on large scales, we have 



'UD,(i^UD,j) 
{t>(1) r (1)' 



TT 



TT 



B 



(i)' 

UD 



C + 2*. 



(101) 
(102) 

(103) 



Thus from Eqs. (|101[) and the continuity of B^jj, ( and '5 across the matching surface, we have 



l (2) ' 
h n p 



,(2)/ 
l ij P 



(104) 



(2) 

The tensor metric perturbation h\- and its time derivative are thus continuous across the matching surface in either 
uniform-density or Poisson gauge for adiabatic perturbations on large scales. 



2. Non-adiabatic matching 

In this subsection, we will consider the non-adiabatic case where the transition surface is determined by a fixed 
value of the field x- Therefore we match on a uniform-x hypersurface which, for non-adiabatic perturbations, need 
not coincide with a uniform-density hypersurface. 

Continuity of the induced metric on the uniform-^ hypersurface implies 

+ = 0. (105) 



,(2) 
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Because we choose E x ^ = 0, Eq. (|87j) implies that B x ^ is continuous across the transition surface. This fact simplifies 
the matching conditions for the extrinsic curvature, Eqs. (1^51) and (j5B"|) , to give 

+ 

(106) 

The second-order gauge transformation for the tensor metric perturbation from a general gauge to the uniform- x 
gauge is given by Eq. (|37p , where from the first-order gauge transformation for a scalar (|2ip gives 





+ 

= -4 




J x_ 







v(l) = 



6 X (1) 



X 



(0) 



(107) 



Thus the matching condition (|105[) can be written in terms of the metric perturbations in an arbitrary gauge as 

+TT 



'h (2) ~ 


+ 

= 2 


[( 









Matching the derivative (|106p in an arbitrary gauge we obtain 



,(2)' 



D 



P)'a«.-fB« - a W) (A (1) 

X.V7 V X / \ ' v 



Ha 



(i) 



+TT 



In the Poisson gauge Bp^ — 0, A^p — ^ and we have 



a ij p 
n ij p 



where 



= -2 



= -4 



MP - 



x x,v 



+TT 



X^J 



+TT 



(1) 



/(0)/ 



H ' 



(108) 
(109) 

(110) 
(111) 

(112) 



Equations (f87|) . (|89|) and (|62|) show that £ x , \l/ and % are continuous, and thus otx^ P is also continuous. Thus Eqs. 
(|110l) and (jllip show that the tensor perturbation, h\- ' p , and its derivative, h\-' , arc continuous in the Poisson 
gauge, even in the non-adiabatic matching case. 

In the uniform-density gauge the situation is rather different. B^ D is given by Eq. {25) and A$ D = (' /H on large 
scales, and we have 



Q 



(1)UD 



X (!) 

OX 



UD 



AO)' 



Cx-C 

H 



Allowing for the fact that ( x and \t are continuous, the matching conditions (|108l) and (|109[) reduce to 



h {2) 

n ij UD 
n ij UD 



= -|[(C + 2*)0 J ]+ TT , 



(113) 

(114) 
(115) 



If £ is not continuous across the transition then the tensor perturbation in the uniform-density gauge, and its time 
derivative, will not be continuous. 

Assuming the evolution is piece-wise adiabatic (that is, £' — before and after the transition) and using Eq. 
for the jump in £ due to the non-adiabatic transition, we have 







+ 2 


'h {2) 








UD_ 


- ~ TP 


h (2) 




+ 4 




UD_ 


- ~ ~H 



4 / w + — W- 



w+ - W- 
l + w+ x ~ 



5 + 3u>. 



-C- \s, 



J x-,n 



H V 1 + w+ 



w+ -w- 1 - 3w_ 

Ov- + ~ ~ <,- > O x - y 



1 + W4 



5 + 3w_ 



PI- 



TT 



(116) 
(117) 
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V. GRAVITATIONAL WAVES FROM INFLATIONARY PERTURBATIONS ON LARGE SCALES 

As an application of the matching conditions for tensor perturbations we will consider the generation of gravitational 
waves on super-Hubble scales after the end of inflation. During slow-roll inflation we have t <1 where the slow-roll 
parameter is defined as e = —H/H 2 = 3(1 + w)/2. 

Using the inflationary solution to set the initial conditions for the tensor metric perturbation on super-Hubble scales 
we can use the solution (|4"u| with the equation of state w_ = — 1 + (2/3)e. During inflation the conformal time ry 
decreases, and we neglect the decaying solutions proportional to Zij p _ and Then the large scale solution in the 
Poisson gauge, Eq. PS)) , is 

h ?h- = y Hp- ~ ^- [diV-djVJT . (118) 

The constant solution on large scales, Yij p _, is the usual free part of the gravitational wave solution whose amplitude 
is determined by the quantum vacuum on small scales, assuming inflation is sufficiently long-lived. We are specifically 
interested in the production of second-order tensor perturbations from first-order scalar perturbations, ( or vp, during 
inflation. Hence we choose to set Yij P _ equal to zero and study the solution 

h ?h- = ~% W-W-r , hf/ p _ = | [aiV-djV-f 7, . (no) 

The tensor metric perturbations produced by the first-order scalar metric perturbations during inflation are thus 
suppressed in the Poisson gauge, and vanish in the limit e — > 0. 

One can use the gauge transformation equation (|51l) during inflation to find the inflationary initial conditions in 
the uniform-density gauge from the Poisson gauge initial conditions as 

h fiuD- =-Jp [diV-d.V-r , hf/ UD _ = 1 (1 - e) [d^d 3 VA TT ■ (120) 

We see that the tensor metric perturbations during inflation produced from first-order metric perturbations are not 
slow- roll suppressed in the uniform-density gauge, in contrast to the Poisson gauge result (|119p . 

Equations (|119j) or (|120[) together with the matching conditions obtained in Section HVl then set the initial conditions 
on large scales at the start of the radiation era that follows the inflationary period. 

A. Adiabatic matching 

First-order scalar metric perturbations can be calculated in any gauge, but the curvature perturbation in the 
uniform-density gauge, remains constant after inflation for adiabatic perturbations and on large scales jB5j, and 
hence we have V+ — V_. The curvature perturbation in the longitudinal gauge, ^, is also continuous for adiabatic 
perturbations, as shown by the matching conditions Eq. (|78|) . but the continuity of ( implies from Eq. (|29l) that "J/' is 
discontinuous at a sudden change in the equation of state. Hence ^ becomes time-dependent after inflation. 

The adiabatic matching conditions for the scalar metric perturbations (|78p determine the amplitude of the decaying 
mode (|8"Tj) in a radiation era, with equation of state w+ = 1/3, following slow-roll inflation 

W + = \ (^-\^%*V.. (121) 

During the radiation era the general solution for tensor metric perturbations in the Poisson gauge, Eq. (|46l) . reduces 
to 

h f P+ = *0>+ + z Hp+V+ 1 + I [d % V+d 3 V + ] TT V 2 + + 2 [cW+cW + ] TT T7+ 4 . (122) 

Using the adiabatic matching conditions for the second-order tensor perturbation in the Poisson gauge (|104l) . which 
require that both the tensor perturbation and its first derivative are continuous at the end of inflation, we find that 
during the radiation era following inflation the integration constants are given by 

Y ijp+ = -4e»#, [d t V-d 3 V-f T , Z ijp+ = -y (l - |e) ^ [diV-djV-f . (123) 
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The general solution for the tensor perturbation in the uniform-density gauge ((50)) for W+ = 1/3 becomes 

h S > UD+ = Y *wd + + ZiiuD + V+ + \ W+diV+T* rf + . (124) 

Either from the adiabatic matching conditions for the tensor perturbations in the uniform-density gauge (|97[) and 
()100p . or from the gauge transformation Eqs. (|S"2"|) of the Poisson gauge results (|123[) . one can obtain the uniform- 
density integration constants as 

Y*vd+ = [diV-diV-f* , Z lJUD+ = -| (l - | c ) r,%. [d i V-d j V-] TT . (125) 

(2) 

Thus the full solution for the induced tensor perturbations, h\j, on large scales during a radiation era after inflation, 
is given in the Poisson gauge (|122[l and uniform-density gauge (|124|) by 

[^y_^y_] TT , (126) 



l v p+ 
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[dtV-djV-] 11 . (127) 

We see that although the time-dependent parts of the tensor perturbation in the two gauges are different, the constant 
mode after inflation in both gauges is the same and it is slow-roll suppressed, Yij p , = Yij UD+ = 0(e), even in the 
uniform-density gauge where the tensor perturbation is not suppressed during slow-roll inflation. 



7 (2) 

y ud+ 



B. Non-adiabatic matching 



In the presence of non-adiabatic perturbations the scalar curvature perturbation Q needs no longer be constant on 
large scales and may be discontinuous, as shown in Eq. (|89p . As a simple model of non-adiabatic perturbations we 
consider the case where the transition from slow-roll inflation to radiation occurs on a uniform-^ hypersurface with 
curvature C, x ^ From the non-adiabatic matching condition for scalar metric perturbations, Eq. (|90p . we obtain 
the amplitude of the constant mode and the decaying mode of ^ in the radiation era following inflation 



w+ 



1 - - 2 ) &- 



5"- 



2eV-)<n: 



(128) 
(129) 



For C x = V- we recover the adiabatic case (| 12 1[) . 

The tensor perturbation in the Poisson gauge, and its time derivative remain continuous for a non-adiabatic tran- 
sition. Matching the inflationary solution ()119|) to the radiation era solution (|122|) with the integration constant V+ 
and Z + given by Eqs. (|128l) and (I129[) determines the two remaining integration constants in the Poisson gauge 



n p+ 



+* 



diV-djV- - 2diC x -djV- 



TT 



Z, 



-877 



+* 



2 in 

ediV-djV- + -(1 - e)8iC x -djC x - - -ed^-djV. 



'FT 



(130) 
(131) 



The tensor metric perturbation in the uniform-density gauge can be discontinuous, with the jump in h^j and its 
time derivative being given by (|116[) and (|117l) . where the incoming non-adiabatic perturbation is given by 



S x - = Cx- - V- • 



(132) 



As a result, matching the inflationary solution (|120[) to the radiation era solution (|124p with the integration constant 
V+ and Z + given by Eqs. (|128[) and (| 129[) gives the two remaining integration constants 



Y, 



i-OUDA 



Jl 3UD+ 



Y, 



-8t?: 



+» 



ediV-djV- 



1 7 

-(1 - ejfliCx-^Cx- ~ Q ed iC x - d J V - 



TT 



(133) 
(134) 
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(2) 

Therefore the full solution for the induced tensor perturbations, h\ ■ , on large scales during a radiation era after 
inflation in this non-adiabatic model for the transition, is given in either the Poisson gauge (I122p and uniform-density 
gauge tHH) by 



i (2) 
v p+ 



n 2 + [(1 - e)d l C x -d J C x - + ed^-djV-]^ + W+* [dtV-diV- - 



TT 



2 1(1 

ediV-djV- + -(1 - 6)0^-8^- - -e^C x -^F_ 



TT 



(135) 



, (2) 
*J [/£>+ 



r) 2 + [(1 - ejaiCx-^iCx- + eftCx- + 4er?2 [fl^K. - 2d i C x -d j V-] 



TT 



n TT 



(136) 



In the adiabatic case where £ x _ = V_ and hence 5 X _ = 0, this reduces to Eqs. (|126[) and (|127[) . The presence of 
non-adiabatic perturbations changes the amplitude of the decaying mode, proportional to Z^j, and the amplitude of 
the growing mode proportional to rj 2 ^. But the free part of the gravitational field, which is gauge- independent, remains 
constant on super-Hubble scales and remains suppressed in the slow-roll limit 



(137) 



VI. CONCLUSION 



Tensor perturbations of the FRW metric become gauge-dependent at second and higher order jl9j . It is however 
possible to construct gauge-invariant tensor perturbations by eliminating the gauge degrees of freedom [24{ . As 
examples we have given gauge-invariant definitions of the tensor perturbation at second order in the Poisson gauge 
and the uniform-density gauge, and the (gauge-invariant) difference between the two. All of this should be familiar 
from the study of gauge-dependent scalar perturbations which are already gauge dependent at first order, and one 
can choose to work in terms of the curvature perturbation in the Poisson (or longitudinal) gauge, 'J, or in the 
uniform-density gauge, £, or both. 

We have focussed on the question of matching conditions for the tensor perturbation. At first order, the Israel 
junction conditions require that the tensor metric perturbation and its time-derivative are always continuous across a 
spacelike hypersurface in any gauge. But at second order we need to consider the behaviour in different gauges. For 
adiabatic perturbations the matching hypersurface should coincide with a uniform-density hypersurface. The Israel 
junction conditions then imply the tensor perturbation in the uniform-density gauge, and its time-derivative, will be 
continuous. For the case of non-adiabatic matching the hypersurface need not coincide with a uniform-density surface 
and so the tensor perturbation in the uniform-density gauge may be discontinuous. However in both cases (adiabatic 
or non-adiabatic) we show that the second-order tensor perturbation in the Poisson gauge, and its time-derivative, 
are continuous. 

As an application we consider the generation of gravitational waves during inflation. We give expressions for the 
general solution for the tensor perturbation on large (super-Hubble) scales in both the Poisson and uniform-density 
gauges when the barotropic index, w is constant. This is a good description of the evolution during the radiation 
dominated era, when w — 1/3, and during a preceding era of slow-roll inflation, when w = — 1 + (2/3)e and e <C 1 
is a slow-roll parameter. It has previously been shown (2|| [27j that the tensor metric perturbation on large-scales 
during inflation is slow-roll suppressed in the Poisson gauge, but not suppressed in the uniform-density gauge. This 
reflects the fact that the scalar metric perturbations are slow-roll suppressed during inflation in the Poisson gauge 
with respect to those in the uniform-density gauge. The question is then, whether the amplitude of gravitational 
waves after inflation is slow-roll suppressed, or indeed whether large tensor perturbation can be generated at second 
order by the transition to a radiation dominated universe. 

Continuity of the tensor perturbation and its first derivative is sufficient to show that the tensor perturbation 
remains slow-roll suppressed after inflation in the Poisson gauge. In the uniform-density gauge tensor perturbations 
are not slow-roll suppressed, reflecting the gauge- dependence of the second-order tensor perturbations. However, the 
constant mode on large scales, Yq, corresponding to the free part of the tensor perturbations is the same in the 
Poisson or uniform-density gauge in the radiation era. It is this constant mode that describes the gravitational waves 
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produced by inflation on large scales. In the case of non-adiabatic perturbations this amplitude can change in the 
uniform-density gauge due to non-adiabatic matching condition, but even in this case its amplitude vanishes in the 
slow-roll limit e — > 0. 

In addition to the constant mode there is a growing mode proportional to rj 2 on large scales in the radiation era. 
This is not the free part of the gravitational field, i.e., it vanishes in the absence of first-order scalar perturbations. 
But, like the scalar metric perturbations during the radiation era, it is not slow- roll suppressed in either Poisson or 
uniform-density gauge. This becomes the dominant term in the tensor perturbations on large scales and corresponds 
to the production of gravitational waves from density perturbations in the radiation dominated era. Ananda et 
al. Q showed that the amplitude of tensor perturbations reaches a maximum at Hubble-crossing, k-q » 1 (where 
our large-scale approximation breaks down), and thereafter follows the usual behaviour for free gravitational waves, 
\h\ oc a -1 , as the scalar metric perturbations oscillate and decay on sub-Hubble scales, krj ^> 1. The first-order 
scalar perturbation, $, decays rapidly on sub-Hubble scales leaving an unambiguous prediction for the amplitude of 
gravitational waves. 

Thus we conclude that the power spectrum of the stochastic gravitational wave background generated from scalar 
metric perturbations from inflation is due primarily to the production of tensor perturbations around Hubble-entry 
during the radiation era, and the second-order gravitational waves from inflation are suppressed in all models of 
slow-roll inflation. 
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Appendix A: Inverse metric up to second order 



From Qfixg = 8" we can obtain the inverse metric order by order. At zeroth order we get 

g°°W = -a- 2 , ,g w (°)=0, = a-*** . 

At first order, it gives 

g 00W = 2a - 2 AW, = a~ 2 B^\ gW> = a~ 2 (2C«<^ - C««) 

The second order components of the inverse metric are 

B (i)^(i)_ 4 ^(i)) 2 + A (2) 



ff oo(2) = -a 



5 o l(2) = a -2 



7 Jy ' — —a 



£»(i)_gj'(i) _ 4 (c^Y S ij + 4C (1) C i:,(1) - C^C^^ - C^S ij + -C ij{2) 



(Al) 
(A2) 

(A3) 
(A4) 
(A5) 



where the indices were raised with 8 lJ . 
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